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1 Introduction 

Let K he a. field of characteristic 0. We denote the differential operator ^ by 
di- The ring of differential operators with polynomial coefficient 

[/3eN" J 
is denoted by and that with formal power series coefficient 

^ = < 5] akA^)d^ I akA^) G K[[x]] I 

[/3GN" J 

is denoted by V. Here x — (xi, • ■ • , Xn) and d — (di, ■ ■ ■ , 

Division theorems are fundamental in several construction for Z?-modules. 
Castro gave a constructive division theorem in 2?, which gives unique quotients 
and remainder (0). However, the division procedure needs infinite reductions. 
On the other hand, Granger, Oaku and Takayama gave a division algorithm for 
algebraic data in T> {^,^), which is an analogous algorithm with the Mora 
division algorithm in the polynomial ring (POj). We call this division the Mora 
division algorithm in D. The division algorithm stops in finite steps, but the 
remainder is not completely reduced. By using the division algorithm, we can 
get a Grobner basis and solve the ideal membership problem for algebraic data 
in V. 

The local b function of a polynomial f G K[x] is the minimum degree monic 
polynomial b{s) G K[s] which satisfies 

3P G V[s] such that Pf+^ = b{s)f'. 

Oaku gave algorithms computing the local b function of a given polynomial by 
using a Grobner basis method in D (jnijEljCZI)- In this paper, we propose 
a new algorithm to compute the local b function. We use the Mora division 
algorithm in D and an approximate division algorithm in P, which gives an 
approximation of the remainder by Castro's division in V. 



2 Division Theorem and Approximate Division 
Algorithm in the Ring of Power Series 

There are several kinds of division theorem in _ftr[[a;]]. Among them, we need 
a division theorem which gives unique quotients and unique remainder for our 
approximate division algorithm in D. We start with explaining this kind of 
division theorem. As to details about the division theorem including a history, 
we refer jTj. 

We define a monomial order <r on monomials in if [[a;]] by the following 
1 X n matrix and a term order < as the tie-breaker. 

Xi • • • 

~Zl 77. ~ 

( term order < ) 

In other words, we define the order <r as 

„ ^ j -ai Q!„ < /3„ or 

I (— ai — — = — /?! — • • • — /3n and x" < a;'^) 

Since the order is not a well order, ordinary division algorithm does not generally 
stop in finite steps. 

Let a(l),--- ,a(s) e (Z>o)". For (a(l), • • ■ ,a{s)), we define a partition 
Ai,--- ,A^A of (Z>o)" as 

Ai =a(l) + (Z>or, 
A2 = (a(2) + (Z>or)\Ai,... 
A^ = ia{s) + (Z>o)") \ (Ai U • • • U A^-^), 
A==(Z>o)"\(Aiu---UA''). 

Let / e ^[[x]]. We define the following notations. 

LM<,, (/) : leading monomial of / with respect to <r 
LT<,, (/) : leading term of / with respect to <r 
LE<^(/) : leading exponent of / with respect to <r 
Exps(/) : set of exponents appearing in / 
Rest<^(/) : /-LT<^(/) 

For example, let / — 3a;i-|-a;ia;2, thenLM<^(/) — xi,LT<^(/) = 32;i,LE<^(/) = 
(l,0),Exps(/) = {(1,0), (1,1)}, andRest<^(/) =xiX2. 

Lemma 2.1. (Division by monomials in if [[a;]]) 

Let a(l),--- ,a(s) S (Z>o)" and / G K[[x]] and A^,--- ,A'',A be the par- 
tition of (Z>o)" with respect to (a(l), • • ■ ,a{s)). Then, there exist quotients 
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qi, - ■ ■ ,qs G -f'^iN] a-iid remainder r G which satisfies the following con- 

ditions : 

a(i) +Exps(g,) C A\ 
Exps(r) C A. 

Especially, (gi, • • • ,qs,r) is uniquely determined by (/, x"^^^ , • • ' i a;"^"-* , <r)- 

Theorem 2.2. (Division theorem in iir[[x]]) 

Let / e gi,-- ■ ,gs G ^[N], and A^, ■ • ■ , A*, A be the partition of 

(Z>o)" with respect to (LE<^(5i), • • • ,LE<^((7s)). Then, there exist quotients 
gi, ■ • • , (7s G ^[[a^]] and remainder r G which satisfies the following con- 

ditions. 

/ = qi9i H — + qsgs + r 

LE<^((7,) +Exps((7,) C A' 
Exps(r) C A 

Especially, (gi, • • • ,qs,r) is uniquely determined by (/, .gi, • • • ,5s, <r)- 
We show a procedure to obtain (?i , • ■ • ,qs,i". 

K [[x]] -division (/, gi, ■ ■ ■ ,gs) 

Input : f eK[[x]],gi,--- ,gseK[[x]] 
Output : (7i, ■ • • ,qs,r G K[[x]] which satisfies 

/ = qigi H — qsgs + r 

LE<„(5,)+Exps((7,) C A* 
Exps(r) C A 

2 : while {F ^Q) { 

3: [g',r'] ^mono-div(F,5i,--- ,5s) 

4 : ^ <7,: + (1 < « < s) 

5 : r <— r + r' 

6: ^^^-E:=i'?^Rest<„(5,) 

7: } 

8 : return [g, r] 

Here, the procedure mono-div(i^, G) at the 3rd line returns q'l,--- ,qs,r' G 
which satisfies 

F = g'iLT<^ (.gi) + • • • + g:LT<„ (.gs) + r' 
LE<„(5,)+Exps(gOc A^ 
Exps(r') C A 

Namely, divide F by LT<^(G) (Lemma EU . 
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This procedure generally does not stop in finite steps. Therefore this proce- 
dure is not an algorithm in a strict sence. If input / and Qi are polynomials, 
we can give an algorithm which returns approximate quotients and approximate 
remainder which are correct up to given total degree N . Approximate algorithm 
of this kind has been discussed by several authors, for example, see ^J. After 
the manner of the approximate algorithm, we will give an algorithm which gives 
approximations of quotients and remainder by ^[[xJJ-division. 

Algorithm 2.3. (iir[[a;]]-approximate-division) 



ii'[[.T]]-approximate-division(/, 51, • • • , gs) 
Input : / e K[x\,gi, ■ ■_,gs G K[x],N E Z>o 

Output : qi, ■ ■ ■ , ql, r, F (z K[x] which satisfies the following conditions. 

^= 9131 H ^qsgs + r + F 

F^0=> |LE<,,(F)| > N 

Terms of g7 whose total degree is smaller than N — |LE<^(5i)| agree 
with those of qi. 

Terms of r whose total degree is smaller than N agree with those of r. 
Here \a\ is J27=i '^i^ where a — (ai, • • • , a„), 
and qi,r are quotients and remainder of _fir[[a:]]-division(/, G). 
T': F ^ f,Ji ^ (1 < i < s)^r ^ 
2 : while (F 7^ and |LE<,,(F)| < iV) { 
3: [g',r'] ^ mono-div(F,G) 

4: qi + q'i {'^ <i < s) 

5 : r <— r + r' 

6: F^-E.=i'Z,'Rest<^(g,) 
7: } 

8 : return [q^r,F] 

(Proof) We denote variables F,'gi,r,q'^,r' at the end of the execution of 
the fc-th while- loop by F^^\'gi^''\r^''\ q'^''\r'^''\ From the result of mono- 
div(F^''~^\G), it holds that 

F''-''=±q'}'^g. + r'^'^+F^'' (1) 

i=l 

LM<,(Qf'^5.) <r LM<^(F^'''^),LM<„ LM<^(f''"'^) (2) 
We will prove that the algorithm stops in finite steps. From ^ and we 

get 

LM<^(f'''') >r LM<,.(F^'"^'^) >r LM<,(F^*-'+'^) >.•••. 
From the definition of the monomial order <r, we get 

|LE<,XF^''')| < |LE<^(F^^+'')| < |LE<^(F^''+'^)| < ■•■ . 
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Since the number of monomials which has the same total degree is finite, there 
exists k such that |LE<^(f )| > Af. This is the condition which stops the 
while loop. Therefore the algorithm stops in finite steps. 

Next, we will prove that the returned qi,r satisfy the conditions of output. 
At the end of the execution of the fc-th while-loop, it holds that 

Since q'^''^ , r'^''^ is the result of mono- div(F^*' ^\g), it holds that 

LE<^(5,)+Exps((zf)) C A' 
Exps(r'('=)) C A 

From these properties, we have 

LM<,(F^'-'^) = max(LM<,(gf • • • ,LM<,(gf )5,),LM<,(r'«)) 

That is 

So we get 

|LE<^(qf ')| > \LE^AF^'-'^)\ - |LE<,(50I, 
|LE<,(r'«)| > |LE<,(F^'-'^)|. 

Since gT^*^' = gl^^^-^' + q'^''^ and r^^^ = r^^^^' + r'C^^ terms of gT^^) whose 

(/c — 1) 

total degree is smaller than |LE<^(i^ )| — |LE<^(5i)| does not change, and 

terms of f^*^^ whose total degree is smaller than |LE<^ (F^'' ) | does not change. 
Therefore qi, f satisfy the condition. ■ 

Example 2.4. (Example of ii'[[a;]]-approximatc-division) 

The example of the case of / = x,G = l + x,N = 5 (execute 
approximate-division(a:, 1 + x, 5)) 



mono-div 




q' 


r' 


Q 


r 


F 
















X 


x = x-l + 




X 





X 





—X ■ X 


-x^ = -x^ ■ 1 + 





-x^ 





o 

X — X 





-{-X^) ■ X 


= • 1 -1- 




x^ 





X - x'^ +X^ 





—X^ ■ X 


-a;4 = -a;4 • 1 + 





-x^ 





X — X + X — X 





— (— X^) • X 


a;5 = • 1 + 




x' 





X — x"^ + x'^ — x"^ + x^ 





—X^ ■ X 



The output of the approximate division algorithm is 

X = {x - x'^ + x^ - x'^ + x^) ■ {1 + x) - x^. 
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In this case, the quotients and the remainder by i^[[x]]-division are 

X = • (1 + .t) + 

1 + X 

= {x-x^ +x^ -x^ ^x'^ ) • (1 + x) + 

3 Division Theorem and Approximate Division 
Algorithm in T>[y\ 

Castro gave a constructive division procedure in T) His division procedure 

returns unique quotients and remainder. In this section, we give an algorithm 
which gives approximations of the quotients and the remainder by his division. 
The approximate quotients and remainder are correct up to a given total degree. 

Let 2/ be a parameter, and ? be a commutative variable which stands for 
d. We define a monomial order < on T>[y\ by using the following 2 x (2n + 1) 
matrix and the tie-breaker <i. 

Xl ■■■ Xn y £,1 ■■■ in 

••• 1 1 ••• 1 

-1 ••• -1 ••• 

( term order <i ) 

And we define an other monomial order <r on V[y\ by using the following 
1 X (2n + 1) matrix and tie-breaker <i. 

Xl ■ ■ ■ Xn y 6 • ■ • in 

-1 ■•• -1 -1 -1 ••• -1 

( term order <i ) 

We define a weight vector e for monomials {x'^^^y^} as 

Xl ■■ ■ Xn y 6 ■ ■ • in 

( ••• 1 1 ••• 1 ) 

If |/?| -|- I7I = \(3'\ + \^'\ (in other words, monomials have the same order with 
respect to the weight vector e), then it holds that 

x"i^y^ < x'^'fy'^' ^ x'^i'^y'' <r x'^' i^' y^' 

Therefore, we get 

LM<(P) = LM<^(ine(P)) 

Here ine(P) is the initial part of P with respect to the weight vector e. For 
example, when P = xidf + x^df + Xidi + 1, we have ine(P) = Xiif + xlif. And 
LM<Jine(P)) = xiilLM^iP) = xiif 

We review Castro's division procedure in Pfy] . Of course, since the monomial 
order < is not a well order, the procedure needs infinite reductions. 
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Theorem 3.1. (Division theorem in Vly], [S]) 

Let P, Pi, • ■ • , P, e V[y], and A\ • ■ • , A", A be the partition of with 

respect to (LE<(Pi),--- ,LE<(Ps)). There exist Qi,-" e 'D[y] which 

satisfies the fohowing conditions. 

P = QiPi + ■■■ + QsPs + R 

LE<(P,) + Exps(g,) C A' 
Exps(P) c A 

Especially, (Qi, • • • ,Qs,R) are uniquely determined by (P, Pi, • • • , Ps, <). 
We show procedure for the division theorem. 

P-division(P,Pi,--- , Ps) 

Input : P,Pi,-- - ,P. ePM 

Output : Qi, • ■ • , Qs, P e I?[y] which satisfies 



P = QiA + • • ■ + OsPs + R 
LE<(P,)+Exps(Q,) C A' 
Exps(P) c A 

1 : Q,; ^ (1 < z < s),P ^ P 

2 : mo <— orde(P) 

3: for (fc ^ mo;fc > 0;fc ^ fc - 1) { 

4 : r <— (total symbol of the part of R whose e-order is k) 

5: [g-,r'] ^ if [[x]] -division (r, ine(Pi), • • • ,ine(Ps),<r) 

6 : Qi ^ (replace ^ with 9 in g-) (1 < « < s) 

7 : P ^ P - EU Q^P^ 

8 : ^ + Q', il<i<s) 
9: } 

10 : return [Qi,R] 



Here, orde(P) is e-order of P, defined by 

orde(P) = max{e • a | a G Exps(P)} , 

and the total symbol of P{x,y,d) = Ea'^a{x,y)d" {aa{x,y) S iir[[x]][?/]) is 
P(a;, y,C) = ^a(a;, For example, when P = (a;i x 1X2)81 + Xidi -I- 1, 

orde(P) = 2 and the total symbol P(x,f) = (xi + 2:ia;2)^i + xiS^i + 1. 

Since the procedure uses if[[x]]-division at the 5th line, it does not stop in 
finite steps. We suppose that inputs are algebraic data, in other words, inputs 
are elements in D. By replacing the i<r[[a;]]-division with _fC[[a;]]-approximate- 
division f Algorithm I2.3|l . we can get approximations of the quotients and the 
remainder by P-division. 

Algorithm 3.2. (2?-approximate-division) 
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X'-approximatc-division(P, Pi, ■ ■ ■ , P,, A'^) 
Input : P,_Pi, • • • ,_f^ (^D[y],N G Z>o 
Output : Qi, • • • , Qs, -R G satisfies 

Terms of Qi whose total degree is smaller than N — |LE<(Pi)| agree with those of Qi 

Terms of R whose total degree is smaller than N agree with those of R 
Here, Qi and R arc the result of X'-division(P, Pi, • • • , Ps). 



1 : Q^ ^ (1 £ z < s),R^ P 

2 : Too <— orde(P) 

3: for (fc ^ 0;fc < TOo;A: <- fc+ 1) 

4 : Mfe ^ max(|LE<(Pj)| + 2(fc - orde(Pii)) | i such that k > orde(Pi)) 

5 : Bound + J2t°o 

6 : for {k ^ toq; fc>0;fc<— fc — 1) { 

7 : f ^ (total symbol of the part of R whose e-order is k and whose 

total degree is smaller than Bound) 

8 : [q^, r'] ^ if [[a;]] - approximate-division (f, ine(Pi), ■ • • , ine(Ps), <r, 

Bound) 

9 : Qi ^ (replace ^ with d in q'^) {1 < i < s) 

10: ^^^ELoOl^^ 

11 : Qi^Qi + Q'i (1 < i < s) 

12 : Bound Bound -Mfc 
13: } __ 

14 : return [Qi,R] 



Since ^^[[ajjj-approximate-division stops in finite steps, the algorithm stops 
in finite steps. In the remainder of this section, we will prove the correctness of 
the algorithm. 

We put 

D{m) = {P& D[y] (or V[y]) \ orde(P) = to} , 

T{i) = [P^ D[y] (or V[y],K[x, (,],K[[x]][y, (,]) \ (the total degree of every term of P) > i 

It follows from the definition that P — Q ^ T{i) is equivalent to that P and 
Q are the same up to total degree i — 1. From the Leibnitz rule, we get the 
following property of the multiplication of approximate elements. 

Lemma 3.3. (Multiplication of approximate elements) 

(T(i) n D{m)) ■ T{j) c T{i + j - 2m) 

Here, dot • means multiplication in the ring of differential operators. 

(Proof) We suppose that P G T{i) n D{m) and Q G T{j). Then we have 

|/3|<m,|a| + |/3|>i |7| + I'5|>J 
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To compute the multiplication PQ, we use the Leibnitz rule. We can get the 
following total symbol for PQ. 

1^ ' ^ (a,/3)GExps(P) (7,5)GExps(Q) 
(Q,/3)eExps(P),(7,/3)GExps(Q),iy<;3,!y<7 

Here, v < (3 means Vi < (3i for all i. 

We consider the minimum total degree of terms of PQ. 

min(|a| + I7I - + |/3| + |<5| - \ {a, P) G Exps(P), (7, 5) G Exps(g), y<P,i^<i) 

(3) 

Since |a| + \j3\ > i, \j\ + \S\ > j and < < orde(P) — m hold, the minimmn 
total degree Q is more than or equal to i + j — 2m. Therefore, we conclude 
PQ eT{i + j ~2m). m 

Lemma 3.4. Let qi and r be the quotients and remainder of 
i^[[x]]-division{/, {51, • • • , <?4 , <,). If / G TiN), then G T(iV-|LE<^(gO|), r G 
T{N). 

(Proof) It holds that 

/ = 9151 H ^ ^sffs + r 

LE<^((7,) +Exps((?0 C A' 
Exps(r) C A 

, where A^, • • • , A*, A are the partition oil^^^ with respect to (LE<^((7i), • • • , LE<^(gs 
From these properties, we get 

LE<^fe5^) <rLE<^(/) 
LE<^(r) <r LE<^(/) 

From the definition of the monomial order <r, we have 

|LE<^(g,)| > - |LE<^(<7,)I ^ q^ e T{N ~ |LE<^(5,)|) 

and 

|LE<^(r)| > |LE<^(/)| = TV ^ r G T{N). 

■ 

Lemma 3.5. Let f,gi,--- ,gs G and / be the part of / whose total 

degree is less than N, that is, f — f E T{N). Let qi and r be the quotients and 
remainder of ii'[[2:]]-division(/, • • • ,gs, <r)- And let ql and r be the quotients 
and remainder of i4r[[x]]-approximate-division(/, gi, • • • ,gs, <r,N). Then qi ~ 
qi€T{N~ |LE<,(g,)l), r-re T{N). 
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(Proof) Let and r' be the quotients and remainder of 
-ftr[[x]]-division(/, gi, • • ■ ,173, Since if [[a;]] -division uniquely gives the quo- 
tient and remainder, the quotients and remainder of if [[x]] -division (/—/, • • • , g 
are qi — q[ and r ~ r' . From f — f E T{N) and Lemma 13.41 qi — q[ G 
T{N - |LE<^(g,)|) and r - r' e T{N) holds. 

From the property of i4r[[a;]]-approximate-division, q'^—qi G r(A^— |LE<^ (.9i)l) 
and r'-r G T{N) hold. So we get q,-q- G r(7V- |LE<^(.gj)|) and r-r G T{N). 

■ 

(Proof of Alaorithm \S. jjj) We will compare steps in the for-loop in P-division 
with those in I?-approximatc division. We put 

Mfe = max(|LE<(Pi)| + 2{k - orde(Fj)) | i such that k > orde(Pj)). 

We suppose that R — R E T{N) in the (fc — l)-th step. Let us compare fc-th 
steps in the for-loop in P-division and 2?-approximate-division. 

2?-division 

r ^ (total symbol of the part of R whose e-order is k) 
[q[,r'] ^ xi[.T]]-division(r, {ine(Fi), • • • , ine(P,)} , <,) 
Q',^qiix,d) 

Rnew ^ R Qi-^i 



2?-approximate-division 

r ^ (total symbol of the part of R whose e-order is k and total degree 
is less than N) 

[q'ijr'] ^ iir[[x]]-approximate-division( r, {ine(Pi), • • ■ ,ine{Ps)} , <r, N) 

We will prove that Rnew— Rnew G T{N — Mk) holds after these computations 
have been done. 

It follows from Lemma 13.51 that 

q[-^ET(N- |LE<^(ine(P.))|) = T{N - |LE<(F,)I), 
r' - F G T{N) 

Remainders r and r are e- homogeneous elements of order fc, and ine(Pi) is 
e-homogeneous element of order orde(Pi). So, if q[,q[ 9^ 0, then q[,q[ are e- 
homogeneous elements of order k — orde(Pi)(> 0). And, if r',r' ^ 0, then 
r',r' are e-homogeneous elements of order k. Therefore, if Q'i,Q[ ^ 0, then 
Q[,Q[ have the same e-order k — orde(Pi) and are the same up to total degree 
N - |LE<(Pi)| - 1. In other words, it holds that 

- Q[ G T{N - |LE<(P,)|) n D{k - ovA,{P,)). 

From Lemma I3. 31 we have 

(T(iV-|LE<(PO|)ni?(fc-orde(PO))-r(0) cT(iV-|LE<(P,)h2(fc-orde(P0)). 
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Therefore it holds that 

Q'iPi - QlPi e T{N - |LE<(Pi)| - 2{k - orde(Pi))) 

Since Mfe = niax(|LE<(Pj)| + 2{k - orde(Pi)) I ^ such that k > orde(Pi)), we 
get 



^0;Pi-^Q^P, GT(7V-Mfe). 



Therefore Rnew - Rnew G T{N - Mk). 

The accuracy of approximation decreases by Mk at each step of the for loop. 
To keep the accuracy, we beforehand add J2^o 

■ 

Example 3.6. (Example of P-division and P-approximate-division) 

Let n = 1. We put P = 9^, Pi = {l + x)d+x. We compare X>-division(P, Pi) 

and 2?-approximate-division(P, Pi, 5). 

At first, we show the procedure of X'-division( P, Pi). 

P4- P,mo ^orde(P) = 2 

r <— (part of R whose e-order is 2) = 

[q[,r'] <- if[[x]]-division(r,ine(Pi), <r) 

R^R-Q[Pi = -d-j^ 

r <— (part of R whose e-order is 1) = — ^ 
Wi,r'] ^ iir[[a;]]-division(r,ine(Pi), <r) 
(9i = -T^,r' = 0) 

R <— R — Q'lPi = 1^^ 

r -f— (part of R whose e-order is 0) = -jq^ ~ 
[q[,r'] ^ is'[[a;]]-division(r,ine(Pi), <r) 

Q'l^O 

R^R-Q[Pi = ^ 

The output of 2?-division(P, Pi) is 



= i- d - )((1 + x)d + x) + 



^1 + X 1 + X 1 + x 

The quotient is 

^ d--^ = {l-x + x'^-x^+x'^ ){d-l) 



\+x l+x 
and the remainder is 
-l + x 



= -1 + 2a; - 2a;2 + 2a;3 - 2a;^ + 2x^ 
l + x 
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Next, we show the procedure of 2?-approximate-division (P, Pi, 5). 



P,mo <-orde(P) = 2 
Mo ^ 0,Afi ^ 1,M2 ^ 3 

Bound ^5 + (0 + 1 + 3) = 9 

f <— ( part of R whose e-order is 2 and total degree is less than 9) = 
[gj, r'] <— i<'[[x]]-approximate-division( r, ine(Pi), <r, 9) 

(^=(l-a; + x2 + x^)Cy = -x^^^) 

Q; ^ (l-^+a;2 + x^)d 

R^R- Q[Pi = -x'^d^ - d- x'^d - 1 + x- x^ + x^ - x^ + x^ - x^ 

Bound ^ 9 - 3^= 6 

f <— ( part of R whose e-order is 1 and total degree is less than 6) = — ^ 
[q[, r'] ^ i4r[[x]]-approximate-division( r, ine(Pi), <r, 6) 
{q[ = -1 + X - x'^ + x^ — x^, r' = x^S,) 
Q[< 1 + x- x^ +x^ - x'^ 

R^R- Q[Pi = -x'^d^ + x^d -x'^d-l + 2x- 2x^ + 2x^ - 2x'^ + 2x^ - x^ 
Bound <— 6 — 1 = 5 

r «— ( part of R whose e-order is and total degree is less than 5) 

= -1 + 2a; - 2x2 _^ ^x'^ - 2x^ 
[q'l, r'] ^ [[x]]-approximate-division( r, ine(Pi), <r, 5) 
= O,?' r = -1 + 2x - 2a;2 + 2x^ - 2x'^) 

R^ R^ -x^d^ + x^d - x^d - 1 + 2x - 2x^ + 2x^ - 2x^ + 2x^ - x^ 



The output of X'-approximate-division(P, Pi, 5) is 

=((1 -x + x^ + ■■■ + x^)d +{-l + x~x^ +x^ - x% ■ ((1 + x)d + x)+ 
- x^d^ - x^d + x^a - 1 + 2x - 2x^ + 2x^ - 2x^ + 2x^ - x^ 

The correct part of the remainder is 

-1 + 2a; - 2a;2 + 2x^ - 2a;'*. 

And this part is the same with the part of the remainder of 2?-division 

4 Computation of the Local b Function by the 
Approximate Division Algorithm in D[s] 

We will apply the approximate division algorithm in D (Algorithm l3.2() to obtain 
the local b function of a polynomial. In the sequel, we use the parameter variable 
s instead of y, and K denotes the field of complex numbers. 
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4.1 Definition and Algorithm of the b Function 



Definition 4.1. (Global b function) 

For given / e if [a;], we define the global b function as the monic polynomial 
b{s) of the least degree which satisfies 3P G D[s] s.t. P ■ /"^^ = 6(s)/*. 

Definition 4.2. (Local b function) 

For given / G -f4r[a;], we define the local b function at the origin as the monic 
polynomial b{s) of the least degree which satisfies 3P G 2?[s] s.t. P ■ /^^^ = 

Example 4.3. (6 function) 

Let / = x{x + y + 1) ~ x'^ + xy + X. The global b function of / is (s + 1)^, 
and {d^ + dxdy) ■ = (s + 1)2/" holds. The local b function of / is s + 1, 
and T+t+-y9- ' ^ + 1)/' holds. 

Oaku gave algorithms computing the global b function and the local b func- 
tion for a given polynomial (511> CH]> El)- ^ Grobner basis method in D is 
used in these algorithm. 

Noro gave an efficient algorithm computing the global b functionf|13|'). In 
this algorithm, a Grobner basis method in D, and modular method are used. 
Especially, to eliminate variables, normal forms with respect to a Grobner basis 
are used. 

In this paper, we propose an algorithm computing the local b function by 
utilizing an approximate normal form with respect to a Grobner basis. Our 
algorithm is an analogous algorithm with Noro's algorithm. For this purpose, 
we review a Grobner basis and a normal form in D and define an approximate 
normal form. 

4.2 Grobner Basis and Normal Form in V[s\ 

Although T>[s\ is a transcedental object, for ideals in T>\s\ generated by elements 
in D[s\, we can compute a Grobner basis in finite steps by either of the following 
algorithms. 

• Lazard's method in D (using a homogenization) 

• method using the Mora division algorithm in (|S], [5]) 

Definition 4.4. (Normal form in T>[s]) 

For P G ©[s] and a Grobner basis G in T>[s\, the remainder of ■D-division(P, G) 
fTheorem l3.1|) is uniquely determined. We call the remainder the normal form of 
P by G with respect to the monomial order <, and we denote it by NF(P, G, <). 

We note that a normal form does not have a finite representation in general. 

Lemma 4.5. (Ideal membership) 

Let G be a Grobner basis of an ideal 2 in V[s\. For P G 2?[s], P G X is 
equivalent to NF(P, G, <) = 0. 
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Lemma 4.6. (Sum of normal forms) 

Let G be a Grobner basis in 2?[s]. For P^Q e T>[s\, it holds that 

NF(P + Q, G, <) = NF(F, G, <) + NF(Q, G, <) 
4.3 Algorithm Computing X fl /ir[s] 

For an ideal X in I?[s], we propose an algorithm computing a generator of X D 
K[s]. In this algorithm, we use approximate normal forms. We fix a Grobner 
basis G of T. 

For g{s) ~ ais^ + a;„is'^^ + ■ ■ ■ + ^I, li holds that 

g{s) GX^NF(.g(s),G,<) = 

^ a,NF(s', G, <) + a/„iNF(s'-\ G, <) + ••• + aoNF(l, G, <) = 0. 

We compute NF(s%G, <) beforehand, and we compute the minimum I and 
coefficients ai, ■ ■ ■ ,ao(zK,ai^O which satisfies 

a;NF(s', G, <) + a/_iNF(s'-\ G, <) + ••■ + aoNF(l, G, <) = 0. 

Then the polynomial a/s' + ••• + ao is a generator of the intersection. 

However we cannot generally compute normal forms in I?[s] in finite steps. 
The normal form NF(s',G, <) has infinite terms in general. Therefore we use 
an approximate normal form. 

Definition 4.7. (Approximate normal form) 

Let P € T>[s] and G be a Grobner basis of an ideal T in I?[s]. We define the 
approximate normal form of P by G up to total degree N as the remainder of 2?- 
approximate-division(P, G, iV) (Algorithm 123 . We denote it by NF(P, G, < 
,N). By definition, NF(P, G, <) and NF(P, G, <,7V) are the same up to the 
total degree N — 1. 

We suppose that 

a/NF(s', G, <, N) + az_iNF(s'-\ G, <, TV) + • • • + aoNF(l, G, <,N) = 0. 

Since we use approximate normal forms, it is not always true that o/s' + 
aj-is'"^ + ■ • ■ + flo G I. However we can solve an ideal membership prob- 
lem for algebraic data in I?[s] by utilizing the Mora division algorithm in D[s]. 
Therefore we need to apply the Mora division algorithm to the gotten candidate 
o/s' + • ■ • + flQ in order to check if it is a member of X. 

We will explain an algorithm computing X D K[s] by using these ideas. We 
suppose X n K[s] ^ {0}. We put 

i^jv = {(ao,--- ,a,)eK'+^ \ a,NF(s\ G, <, iV) + • • • + aoNF(l, G, <, iV) = 0} 
= {(ao, • • • , a.) e K'+^ \ a,NF{s\ G, <) + ••• + aoNF(l, G, <) = 0} , 
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where G is a Grobner basis of X and iV G N. It holds that 

The polynomial a.is^ + • • • + ag which satisfies Li ^ {0},Li_i = {0} and 
(ao, • • • ,a,;) £ Lj is a generator of Z H In summary, an algorithm to 

compute the generator is as follows: 

Algorithm 4.8. (Computing the generator of X n K\s\) 

(1) N ^ (a natural number), d ^ 

(2) Find i which satisfies Li^i^jq — {0}, L^jv {0}. (Compute Ld^N, Ld+i,N, ■ ■ ■ ■ 
Note that L.^im = {0} ^ Lq = Li ^ ■ ■ ■ ^ L,^i = {0}.) 

(3) Take an element (a^, • • • , oo) G Li,pf and put g{s) ~ a,;s' + ■ • • + ao. 

(4) Divide g{s) by G by using the Mora division algorithm. If the remainder 
is 0, then g{s) is the generator. If not, then set N ^ N + l,d ^ i and 
goto (2). 

4.4 Algorithm Computing the Local b Function 

For a polynomial / £ K[x], an algorithm of computing the local b bmction of / 
at the origin is as follows: 

(1) Compute a set of generators of Ann^jj^j/*. Denote it by H. 

(2) Set I = V[s\{H U {/}). Compute the generator b{s) of X n K[s\. The 
polynomial b{s) is the local b function. 

In (1), we take a set of generators of Annjjjs]/* as a set of generators of 
Ann^j^j/**. We can compute a set of generators of Ann^ij^j/* by Oaku's algo- 
rithm (jni)- In (2), since I n K[s\ ^ {0} holds (from the existence of the local 
b function), we can utilize the Algorithm 14.81 

We show an example of computing the local b function. 

Example 4.9. (Local b function oi ] — x^{y + l)^z^ at the origin) 

The global b function of / is (s + l)^(s + ^)'^, and the local b function is 

(,+ l)2(,+ l)2. 

A set of generators H of Ann.gj^j/'' is 

{Pi = -2s + zd,_,P2 = -xd,j, + P3 = -dy + - ydy\ 
A Grobner basis G di J — T>\s\(Il U {/}) with respect to < is 

{f,Pi,P2,P3,P4 = -xz^d, - 2xz^, P5 = -z^a^ - 2xz^d^ - 2z3c), - 2z^} 
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We set TV = 1 and will apply the Algorithm 14.81 to compute the intersection 
J'r\K[s]. We note that approximate normal forms of s* by G up to total degree 
iV = 7 are 



NF(1,G, <,7) = 1 
NF(s,G,<,7) = ^zd, 

NF{s\G,<,7) = lz'd! + lz'd! + lzd. 
NF(.^G,<,7) = -^z^d! - gz^af - - i. 

Steps of Algorithm 14.81 are 

Li,i = {0} Li = {0} 

^ /n 1 \ Mora division 

L2.1 3 (U, Ij 5 > non-zero 

G 

■J- 1 \ Mora division 

^2,2 3 (u, Ij s > non-zero 

^2,3 - {0} £2 = {0} 

1 1 \ 9 1 Mora division 

^3,3 3 (0, — 1) - -^s > non-zero 

G 

f /r^ 1 1 \ 9 1 Mora division 

^3^4 3 (0, —"1, 1) — > non-zero 

i^3,5 - {0} L3 = {0} 

r /^1 "^-.x '^9 1 Mora division 

^4,5 3 (0, ^, 1) s'^ - ^s^ + > non-zero 

T 3-i\ 3 39,1 Mora division 

£4,6 9 (0, 5, 1) s - fs-^ + > non-zero 

Uj = {0} L4 = {0} 

L5,7 3 (i, I, f , 3, 1) s^ + + lls^ + ls + \ ^l^-^li™ 0. 
Therefore, the local b function is + Ss^ + + |s + ^ = (s + 1)^(5 + i) 



5 Implementation and Timing Data 

Our algorithm computing the local h function has been implemented by utiliz- 
ing computer algebra system "Risa/Asir" ('|19|'). The name of our package is 
"nk_mora/local-b.rr". We show the timing data. 

Our algorithm computing the local h function consists of the following 3 
parts. 

(1) Computation of a set of generators of Ann^j^j/" (In fact, computation of 
a set of generators of Ann^ij^]/'') 

(2) Computation of the Grobner basis of X = Anngr^,/'' + V[s\f 
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(3) Computation of the intersection X n iir[s] f Algorithm 14. 8|l 

To perform the step (1), we use Oaku's algorithm computing Annn[,] f ^f|15)'). 
In the timing table, we denote this by (Z). To perform the step (2), we have 
the following 2 algorithms. 

(2a) Buchberger algorithm utlizing the Mora division algorithm in D[s] 
(2b) hazard's homogenized method in D[s] 

We denote this parts by (GB-a) and (GB-b). To perform the step (3), we use 
the Algorithm 14. 81 We denote this part by (localb-nf). 
We took some examples from 



/ 


I 


GB-a 


GB-b 


localb-nf 


localb 


dog 


X + y'' + z'' 


0.00504 


0.0190 


0.00494 


0.00503 


0.00612 


1 


2 1 2 1 2 

X +y + z 


0.00524 


0.00157 


0.00224 


0.0229 


0.0119 


2 


3 1 2 1 2 

X +y + z 


0.00637 


0.0298 


0.0694 -1- 0.0192 


0.0417-1-0.0244 


0.0226 


3 


x^ + y^ + z" 


0.00631 


0.0317 


1.41 -f 0.27 


0.131 + 0.0884 


0.0384 


4 


a;'^ + + z^ 


0.00624 


0.0306 


7.23 -1- 1.51 


0.325 + 0.272 


0.0612 


5 


a;" + + z^ 


0.00619 


0.0291 


33.0 -1- 3.25 


0.726 + 0.598 


0.129 


6 


a;'' + y2 + z^ 


0.00610 


0.0281 + 0.0170 


137.5 + 24.1 


1.26 + 0.674 


0.247 


7 


x-^ + xy'' + z^ 


0.0260 


0.172 


0.107 


0.114 + 0.0756 


0.0767 


4 


x^ + xy'^ + z^ 


0.0566 


0.140 


2.29 -1- 0.338 


0.501 -1- 0.284 


0.184 


6 


a;^ + xy'^ + z^ 


0.0903 


0.194 


32.6 -1- 5.47 


0.505 -1- 0.284 


0.173 


6 


a;^ + + z- 


0.00858 


0.118 


1.45 


1.43 + 0.821 


0.197 


7 


3 1 3,2 

a;^ + xy' + z 


0.0110 


0.236 + 0.337 


1.09 


2.94 + 1.40 


1.06 


8 


a;^ + y^ + z^ 


0.00979 


0.123 


18.3 + 2.82 


4.36 -1- 1.80 


0.615 


9 


a;" + a^-^y^ + y^ 


0.70 




0.052 




7.06 


7 


a;'' + :E=y^ + y" 


0.0061 


0.0138 


0.0019 


0.163 


0.153 


6 


a;^ + K^y^ + y^ 


0.031 


0.0748 


0.018 




0.0932 


4 


a;^ + y^ + z^ 


0.0062 


0.00162 


0.0025 


0.142 


0.112 


5 


a;^+y" + z^ 


0.0134 


0.052 






23.2 


18 


3 1 2 2 

a; + y z 


0.029 


0.0462 


0.016 


0.868 -f 0.108 


0.144 


7 


/ 3 2 2\2 

(x -y Z ) 


0.102 


7.49 + 1.67 


0.0932 


17.4 + 1.40 


3.99 


14 


x^ - y^z^ 


0.013 


0.0827 


0.092 


3.77 + 0.279 


0.817 


13 


x^ — y^z^ 


0.0088 


0.0265 


0.0028 


15.9 -1- 1.22 


3.53 


17 


x^ + y'^ — 3xyz 


0.019 


0.473 


0.020 


0.641 + 0.0540 


0.0946 


5 


x^ + y^ + z^ — 3xyz 


0.013 


0.140 


0.012 


0.0658 


0.0207 


3 


/ 5 2 2 \ 

y{x -y z ) 


0.014 


16.7 + 3.43 


0.58 


58.6 -1- 5.24 


6.10 


18 


y(x? -y'z') 


0.644 


2.65 -1- 0.668 


0.15 


4.18 + 0.338 


0.711 


10 


y((y+l)a:^-y2z2) 


0.284 


18.8 + 1.34 


3.6 + 0.14 




10.2 
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• localb — total time of Oaku's algorithm computing the local 6 function 

m) 

• deg — degree of the local h function 

• machine — CPU : Athlon MP 1800-1- (1533MHz) (2 CPU), Memory : 
3GB, OS : FreeBSD 4.8 

In the anterior half, Buchberger algorithm with Mora division (GB-a) is 
faster than Lazard homogenize method (GB-b). While, in the last half, (GB-b) 
is faster than (GB-a). From our experiments, in general, it seems that (GB-b) 
is faster than (GB-a). 

In almost every case, Oaku's algorithm(localb) is faster than our algorithm(T 
and GB-* and localb-nf). The reason is that the computation of Grobner basis 
and division in T) is heavy. Especially, the computation of I?-division is heavy. 
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